Quasi-long range order in the random anisotropy Heisenberg model 
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The large distance behaviors of the random field and random anisotropy Heisenberg models are 
studied with the functional renormalization group in 4 — e dimensions. The random anisotropy model 
' is found to have a phase with the infinite correlation radius at low temperatures and weak disorder. 

, The correlation function of the magnetization obeys a power law (m(ri)m(r2)) ~ |rj — r2|~°' 62e . 

The magnetic susceptibility diverges at low fields as \ ~ jy _1 + 015<! . j n the random field model the 
correlation radius is found to be finite at the arbitrarily weak disorder. 

Q 

\Q ] The effect of impurities on the order in condensed matter is interesting since the disorder is almost inevitably 
t-H ■ present in any system. If the disorder is weak, the short range order is the same as in the pure system. However, 
the large distance behavior can be strongly modified by the arbitrarily weak disorder. This happens in the systems 
of continuous symmetry in presence of the random symmetry breaking field |jj . The first experimental example of 
this kind is the amorphous magnet H|j|. During the last decade a lot of other related objects were found. These 
are liquid crystals in the porous media [Q], nematic elastomers He-3 in aerogel |j) and vortex phases of impure 
superconductors [Q. The nature of the low-temperature phases of these systems is still unclear. The only reliable 
statement is that a long range order is absent [f[|]^-|i"(i|] . However, other details of the large distance behavior are 
poorly understood. 

The neutron scattering JlT[ | reveals sharp Bragg peaks in impure superconductors at low temperatures and weak 
I | external magnetic fields. Since the vortices can not form a regular lattice ||, it is tempting to assume, that there is 
a quasi-long range order (QLRO), that is the correlation radius is infinite and correlation functions depend on the 
distance slow. Recent theoretical |Q and numerical studies of the random field XY model, which is the simplest 
model of the vortex system in the impure superconductor Q , support this picture. The theoretical advances jl2j are 
afforded by two new technical approaches: the functional renormalization group [JL4| and the replica variational method 
JL5| . These methods are free from drawbacks of the standard renormalization group and give reasonable results. 
The variational method regards a possibility of spontaneous replica symmetry breaking and treats the fluctuations 
approximately. On the other hand, the functional renormalization group provides a subtle analysis of the fluctuations 
about the replica symmetrical ground state. Surprisingly, the methods suggest close and sometimes even the same 
■ results. 

Both techniques were originally suggested for the random manifolds [jb|Jl5| and then allowed to obtain information 
about some other disordered systems with the abelian symmetry fl^ , |l6| -|l8|jT It is less known about the non-abelian 
systems. The simplest of them are the random field 0] and random anisotropy Heisenberg models. The latter was 
introduced as a model of the amorphous magnet [^) . In spite of a long discussion, the question about QLRO in these 
models is still open. There is an experimental evidence in favor of no QLRO jl9| . On the other hand, recent numerical 
simulations [^O) support the possibility of QLRO in these systems. The only theoretical approach, developed up to 
1 . | now, is based on the spherical approximation pi] . However, there is no reason for this approximation to be valid. 

In this letter we study the random field and random anisotropy Heisenberg models in 4 — e dimensions with the 
functional renormalization group. The large distance behaviors of the systems are found to be quite different. While 
O ■ in the random field model the correlation radius is always finite, the random anisotropy Heisenberg model has a 
phase with QLRO. In this phase the correlation function of the magnetization obeys a power law and the magnetic 
• *h . susceptibility diverges at low fields. 

To describe the large distance behavior at low temperatures we use the classical nonlinear er-model with the 
£3 . Hamiltonian 

_ _ i 

H= /"d D x[J^^n(x)^n(x) + ^ mp (x)], (1) 
J ii 

where n(x) is the unit vector of the magnetization, Vi mp (x) the random potential. In the random field case it has the 
form 

Vimp = -/] fea(x)n ffi (x); a = x, y, z, (2) 
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where the random field h(x) has a Gaussian distribution and (/i a (x)/i J g(x')) = A 2 5(x — x')<$ Q| g. In the random 
anisotropy case the random potential is given by the equation 

V imp = r Q/3 (x)n Q (x)n /3 (x); a,0 = x, y, z, (3) 

a,/3 

where r Q ^(x) is a Gaussian random variable, (t q ^(x)t 7< 5(x')} = A 2 S a7 5psS(x — x'). Random potential (^) corresponds 
to the same symmetry as a more conventional choice VJ mp = — (hn) 2 , but is more convenient for the further discussion. 

The Imry-Ma argument [QJ|] suggests that in our problem the long range order is absent at any dimension D < 4. 
One can estimate the Larkin length, up to which there are strong ferromagnetic correlations, with the following 
qualitative renormalization group (RG) approach. Let one remove the fast modes and rewrite the Hamiltonian in 
terms of the block spins, corresponding to the scale L = ba, where a is the ultraviolet cut-off. Then let one make 
rescaling so that the Hamiltonian would restore its initial form with new constants A(L), J{L). Dimensional analysis 
provides estimations 

J{L) ~ b D - 2 J(a);A(L) ~ b D/2 A(a) (4) 

To estimate the typical angle </> between neighbor block spins, one notes that the effective field, acting on each spin, 
has two contributions: the exchange contribution and the random one. The exchange contribution of order J(L) is 
oriented along the local average direction of the magnetization. The random contribution of order A(L) may have 
any direction. This allows one to write at low temperatures that <j>{L) ~ A(L)/J(L). The Larkin length corresponds 
to the condition <j>(L) ~ 1 and equals L ~ ( J/A) 2 /( 4 ~ D "> in agreement with the Imry-Ma argument 0. If Eq. were 
exact, the Larkin length could be interpreted as the correlation radius. However, there are two sources of corrections 
to Eq. (0). Both of them are relevant already at the derivation of the RG equation for the pure system in 2 + e 
dimensions p^] . The first source is the renormalization due to the interaction and the second one results from the 
rescaling of the magnetization, which is necessary to ensure the fixed length condition n = 1. The leading corrections 
to Eq. ffl) are proportional to 4> 2 J,(j) 2 A. Thus, the RG equation for the combination (A(L) / J (L)) 2 is the following 



d (A{L)V_ (A{L)\ 2 fA(L)\ 4 
d\nL\j{L)) \J{L)J + \J{L)) 



, e = 4-Z) (5) 



If the constant c in Eq. ([3|) is positive, the Larkin length is the correlation radius indeed. But if c < 0, the RG 
equation has a fixed point, corresponding to the phase with the infinite correlation radius. As it is seen below, both 
situations are possible, depending on the system. 

To derive the RG equations in a systematic way we use the method, suggested by Polyakov [[22J for the pure system. 
The same consideration as in the XY jjjj and random manifold |Q models suggests that near a zero-temperature 
fixed point in 4 — e dimensions there is an infinite set of relevant operators. After replica averaging, the relevant part 
of the effective replica Hamiltonian can be represented in the following form 

H R = J d D x[Y^ i £ - £ (6) 

a )i ab 

where a, b are replica indices, R(z) is some function, T the temperature. In the random anisotropy case the function 
R(z) is even due to the symmetry with respect to changing the sign of the magnetization. 

The one-loop RG equations in 4 — e dimensions are obtained by a straightforward combination of the methods of 
Refs. jlij and p2| . The equations below are given for the arbitrary number N of the components of the magnetization. 
In the Heisenberg model N = 3. The RG equations become simpler after substitution for the argument of the function 
R(z): z = cosc/>. In terms of this new variable one has to find even periodic solutions R(4>). The period is 2ir in the 
random field case and 7r in the random anisotropy case. In a zero-temperature fixed point the one-loop equations are 

T =-(D-2)-2(N-2)R"(0) + O(R 2 ,T); (7) 



d\nL 



= ^§ = eJW) + [R"m 2 - 2R"(4>)R"(0) 
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(N - 2)[AR{4>)R"(Q) + 2ctg0i?'(0) J R"(O) - (f^) ] + 0(R 3 ,T) (8) 
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The two-spin correlation function is given by the expression (n(x)n(x')) ~ |x — x' v , where 

r) = -2(N - 1)R"(Q) (9) 



The same equations (|7[j9j) were derived by a different method in Ref. J23|. In that paper the critical behavior in 4+ |e| 
dimensions was studied by considering analytical fixed point solutions R(4>). In the Heisenberg model, analytical 
solutions are absent and they are unphysical for N ^ 3 p3|| . In this letter we search for non-analytical R((f>)- 

As shown below, the random field model can be completely studied by analytical means. In the random anisotropy 
case, one has to solve Eq. (||) numerically. Since coefficients of Eq. (||) are large as <ft — * 0, it is convenient to 
use the expansion of R(4>) over \<f>\ at small (f>. At larger <fi the equation is integrated by the Runge-Kutta method. 
The solutions to be found have zero derivatives at <j> — 0,ir/2. At N — 3 the solution with largest i?"(0)|, which 
corresponds to rj = 0.62e (||), has two zeroes in the interval [0,7r]. There are also solutions with 4 and more zeroes. 
They all correspond to rj < 0.5e. We shall see below, that these solutions are unstable. 

To test the stability of the solution with two zeroes, we use an approximate method. The instability to the constant 
shift of the function R(4>) has no interest for us, since the constant shifts do not change the correlators jF4j. To study 
the stability to the other perturbations, it is convenient to rewrite Eq. (||), substituting uj(R"((j))) 2 for (R"((f>)) 2 . The 
case of interest is u> — 1 , but at u> = the equation can be solved exactly. The solution at u> = 1 can be found with 
the perturbation theory over w. The exact solution at lo = is iC=o(0) = e(cos</>/24 + 1/120). The perturbative 
expansion provides the following asymptotic series for r\: rj = e(0.67 — 0.080-1 + 0.14w 2 — ...). The resulting estimation 
7] = e(0.67 ± 0.08) agrees with the numerical result well. This allows us to expect that the stability analysis of the 
solution i? w= o of the equation with ui = provides information about the stability of the solution of Eq. (|8|). A 
simple calculation shows that R u= o is stable in the linear approximation. Thus, there is a stable zero-temperature 
fixed point of the RG equations with the critical exponent of the correlation function 

r) = 0.62e (10) 

The critical exponent 7 of the magnetic susceptibility x(H) ~ H^ 1 in the weak uniform field H is given by the 
equation 

7 = 1 + (N- l)R"{Q)/2 = l-0.15e (II) 

Let us demonstrate the absence of physically acceptable fixed points in the random field case. We derive some 
inequality for critical exponents. Then we show that the inequality has no solutions. We use the rigorous inequality 
for the connected and disconnected correlation functions [E4| 



(n a (q)n a (-q)) - (n a (q)n b (-q)) < const v / (n Q (q)n a (-q)), (12) 

where n(q) is a Fourier-component of the magnetization, a, b are replica indices. In a fixed point, Eq. (|l^) provides an 
inequality for the critical exponents of the connected and disconnected correlation functions [^ij . The large distance 
behavior of the connected correlation function in a zero-temperature fixed point can be derived from the expression 
X ~ / ^(( n (0)n(x)))d £) x and the critical exponent of the susceptibility (|Tl|). Finally, one obtains the following relation 



N 

±-D<— v , (13) 

where r\ is given by Eq. (g). This equation does not have solutions at N = 3. At N > 3 Eq. (|l^) is incompatible 
with the requirement r\ > 0. Thus, there are no accessible fixed points for N > 3. 

In the previous paragraph, inequality (|l3| ) is derived for the model (^) with the Gaussian random field (0). It can 
also be extended to a more general situation (||) . If one adds a weak Gaussian random field (^J) to any Hamiltonian, 
it suffices for Eq. (|l^) to become valid. The addition of the Gaussian random field corresponds to the transformation 
i?(n a rib) — » i?(n a n(,) + An a n/, in Eq. (||), where A > is a constant. Thus, if at some A the function i?(n a n/,) = 
i?(n a n(,) — An a nfc is possible as a disorder-induced term in Eq. then Eq. ( |l2] ) is valid for the system with the 
disorder-induced term i?(n a nfc). Finally, we conclude, that inequality (|l^) may be broken only for Hamiltonians (|^), 
which lie outside the physically acceptable region or on its border. This suggests the strong coupling regime with a 
presumably finite correlation radius. 

In the random anisotropy case a similar consideration uses the connected and disconnected correlation functions 
of the field (n x (r)n y (r)) in presence of Gaussian disorder (0). The resulting condition for the critical exponent, 
V > {N — l)e/4, rules out all but one fixed points of RG equation (N). 
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The question of the large distance behavior of the random field and random anisotropy Heisenberg models was 
discussed by Aharony and Pytte on the basis of an approximate equation of state pf| . They also obtained QLRO in 
the random anisotropy case and its absence in the random field model. However, we believe, that this is an occasional 
coincidence, since the equation of state |2q ] is valid only in the first order in the strength of the disorder, while higher 
orders are crucial for critical properties |2(|. In particular, the approach ]^5| incorrectly predicts the absence of QLRO 
in the random field XY model and its presence in the random anisotropy spherical model. It also provides incorrect 
critical exponents in the Heisenberg case. 

The random anisotropy Heisenberg model is relevant for the amorphous magnets In the same time, for their 
large distance behavior the dipole interaction may be important Besides, a weak nonrandom anisotropy is 

inevitably present due to mechanical stresses. 

In conclusion, we have found, that the random anisotropy Heisenberg model has the infinite correlation radius and 
a power dependence of the correlation function of the magnetization on the distance at low temperatures and weak 
disorder in 4 — e dimensions. On the other hand, the correlation radius of the random field Heisenberg model is always 
finite. 

The author is thankful to E. Domany, G. Falkovich, Y. Gefen, S.E. Korshunov, Y.B. Levinson, V.L. Pokrovskiy and 
A.V. Shytov for useful discussions. This work was supported by RFBR grant 96-02-18985 and by grant 96-15-96756 
of the Russian Program of Leading Scientific Schools. 
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